“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1974-03 


Hydraulic ram shock phase effects on fuel cell survivabilit 


Kappel, Leslie George 


Monterey, California. Naval Postgraduate School 
http://ndl.handle.net/10945/17060 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist Ser Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ies) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


HYDRAULIC RAM SHOCK PHASE EFFECTS 
ON FUEL CELL SURVIVABILITY 


Leslie George Kappel] 





NAVAL POSTGRADUATE SCHOOL 


Monterey, California 


HYDRAULIC RAM SHOCK PHASE EFFECTS 


ON FUEL CELL SURVIVABILITY 


by 


Leslie George Kappel 





Thesis Advisor: PL DOW 7 ise 


March 1974 


ante 


Approved {or public rckease; distribution wilincted. 








Hydraulic Ram Shock Phase Effects 


on Fuel Cell Survivability 


by 


- Leslie George Kappel 
Lieutenant, United States Navy 
B.S., Miami University (Ohio), 1966 


Submitted in partial fulfillment of the. 
requirements for the degree of 


MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 
March 1974 





ABSTRACT 


Hydraulic ram concerns the dynamic loading and cata- 
strophic failure of liquid filled fuel tanks impacted by high 
speed projectiles. Hydraulic ram is divided into two phases: 
shock phase and cavity phase. The shock phase was studied. 
Theoretical predictions of shock radius versus time were com- 
pared with experimental results. The theory was found to be 
reasonably accurate for blunt shaped projectiles. 

Included in the theory were predictions of Shecmiardssure 
and pressure profiles behind the shock. Theoretical pres- 
sures are presented graphically but were not validated by 
experiment. 

Velocity decay of the projectiles in the liquid was in- 
vestigated and correlated with theory to provide information 


for further studies of the cavity phase. 
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Area of deformed surface of the projectile 
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decay equation 
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Drag coefficient 
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Initial kinetic energy of the projectile; energy re- 
leased at a point source 


Change in projectile energy due to shock formation 
Constant defined by equation (2.14) 
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Shock acceleration 
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Non-dimensional pressure; equivalent flat plate area 
Non-dimensional pressure after the shock 
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Mass of projectile 

Exponent defined in the equation of state 
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I. INTRODUCTION 


Non-nuclear survivability of aircraft flying in a hos- 
tile arena has been studied throughout the history of avia- 
tion. Recent experiences in Vietnam have demonstrated the 
vulnerability of sophisticated aircraft to both small arms 
ground fire and the more complex surface-to-air defenses. 
Except for the rare canes a direct impact by a surface- 
to-air missile, the damage has been incurred by relatively 
small projectiles. | 

Tactical aircraft fuel tanks or cells have the largest 
Surface area and volume of all vulnerable components. For 
this reason fuel cells have been most susceptible to damage. 
Projectile impact with fuel cells have caused fuel starva- 
tion, fire or explosion as well as extensive structural dam- 
age. The En reTvaDTT Tey of aircraft when a ballistic threat 
causes destruction of fuel filled tanks has been the main 
impetus behind this study. Astronomical costs of combat 
aircraft have made it cost effective to acquire more knowl- 
edge.of the destructive mechanism to fuel cells. MnO ner 
application of this information could be in space flight 
where meteoroid collisions are possible. 

Projectiles which penetrate fluid filled cells cause 
damage many times more severe than that incurred by impact 
With an empty cell. With sufficient projectile energy there 
will be catastrophic failure of cell structural components. 


Hydraulic ram is the term used to describe such destructive 





events. The hydraulic ram phenomenon is usually described 
by two phases: the shock phase and cavity phase. The effects 
of these phases are assumed to be sequential, occurring in 
the order listed above. Upon initial projectile impact, 
energy 1S transferred to the fluid, creating a strong hemi- 
spherical shock centered at the point of impact. The shock 
pressure loading may cause failure of the entrance wall. 
Pressures load the wall for a very short time because of 
geometrical attenuation. Essentially there is an impulsive 
load on the front wall during the shock phase. 

As the projectile travels through the fluid, a cavity is 
formed. This phase has sometimes been called the bubble 
phase because of the bubble-like appearance of the cavity. 
Others have called this phase the drag phase, using the la- 
bel to describe the mechanism which creates the cavity [Ref. 
1]. Oscillation of the cavity volume produces pressure 
pulses that may be sufficiently intense to cause catastroph- 
ic failure of cell components. The intensity of the result- 
ing pressure pulSes is weaker than the shock phase pressure 
pulse but the duration is much longer [Ref. 2]. Ultimate 
destruction of tank components other than the entry wall has 
been generally attributed to the cavity phase. 

Until recently, experimental data from studies of hy- 
draulic ram have not been able to describe the baSic phenom- 
enon. Entire aircraft or wing sections have been destroyed 
' while studying the phenomenon. This method of experimenta- 


tion has become too costly; hence theory explaining the 





phenomenon is necessary to defeat hydraulic ram. Much of 
the theory as it applies to hydraulic ram is still in its 
infancy. The complete mechanism is complex, requiring ac- 
curate and detailed analytical efforts. Correlation of 
theory with experiments must be accomplished simultaneously 
to assess theory validity. 

Research at the Naval Postgraduate School has been de- 
Signed to isolate and observe the individual aspects of hy- 
draulic ram. The ballistic range has been used to study 
energy losses through aluminum plates, shock propagation and 
the cavity phase. Specific details were observed with the 
intent of applying the knowledge gained to the design of a 
survivable fuel cell. 

This study has investigated the shock phase usSing the 
theory of Yurchovitch and Lundstrom [Refs. 1 and 2]. A com- 
puter program was formulated from the mathematical model. 
Experimental data for various projectile energies were com- 
pared with the results of the theory. Yurchovitch's formu- 
lation of shock radius as a function of time has agreed with 
the experimental data. These tests also have shown that 
damage to the tank front wall was not extensive enough for 


Significant fuel leakage at the energy levels tested. 
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II. ANALYTICAL STUDY 


In order to predict the effects of the shock phase of 
hydraulic ram, it is necessary to study how the projectile 
impact energy affects the fluid that is disturbed. It was 
assumed that the impact energy, Ej, was known and was suffi- 
cient to create a strong hemispherical shock. The perturba- 
tion in the fluid can then be described by assuming that the 
shock radius is proportional to time to the 0.8 power, and 
that there is a power law density profile behind the |shodl. 
This formulation was suggested by Yurchovitch [Ref. 2]. 
Yurchovitch found that for a given fluid equation of state 
the energy transferred to the fluid determines the resulting 


flow field. 


A. SHOCK FRONT PRESSURES 

The one-dimensional Rankine-Hugoniot equation for a 
shock moving through undisturbed fluid was used’ to determine 
shock front conditions. To complete the solution an equation 
of state was needed. The Tait Equation for isothermal com- 
pressibility of liquids, 
7 Ea(ex) = 4 ’ (25515) 
Vv n[P+B (s) ] 
WNEHEGewneds a Constant for a specific liquid and B is a func- 
tion of entropy only, was integrated along lines of constant 


pressure and entropy. Integrating equation (2.1) yields 


11 





v(T, v(T,0))% 


= B 
(s 2a = 


als (2.2) 


Assuming the changes across the shock are adiabatic and the 
pressure in the undisturbed fluid is small compared to the 


shock front pressure, B iS given as 


c2 
BO SCme here ex (AP) | a) 


n dp E 


With the assumption that B and n are constant for a given 


fluid the equation of state is 


Pofo, (P yn 


—1) | (2.4) 
n Po 





For water n, a constant, was found to be 7.0 [Ref. 3]. 
Non-dimensionalization of the variables was found to be 


convenient and the following scheme was adopted: 


r 


o= R. (t) 
Ms = Rg(t)/C, 


wap (r,t) 
p(e,M,) = 


ee (2.5) 
R(t) 


(E,Mg) 


Par, 3) 


£(E,M.) 
eRe AGE) 


e(r,t) 
Bie a(te) 
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g(E Mg) 
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With the assumption that P,<<P and e,<<e the Rankine-Hugoniot 


oO 
equations become 


ihe wv, (1-9,) a 
f, = >, (2, 6) 
5, a sf, Ge) 


1 


and the equation of state is 


eae 1) (2.7) 
nM¢ 


Solving the above equations for the parameters at the shock 


the following is obtained: 


1 
es (2.8) 
© 1-4 (yh-1) 
nMs 
di = uz etd) (2.9) 
f, = b, (2.10) 


By assuming values for y,, the equation (2.8) can be solved 
for the shock Mach number, M.. Subsequently, f; and $) ver- 
sus M, can be found. Conditions at the shock front are now 


known while conditions behind the shock front are still un- 


known. 


Pees RE OOURMEF EE LD BPH EIND eT riMe SHOCK 


The density profile behind the shock is assumed to have 


the form 


yp = yeas) (2.11) 


3 











To find q(M.) the total mass enclosed by the shock, 
Ri 2n SWE? 
PoRs2m/pe de , 


is equated to the original mass occupied by the enclosed 
volume, 

3 

2/3TRSP, 

with the result 


fyetae — 7 (2.12) 


Substituting the assumed density profile into (2.12), q (M.) 
is found to be 
q(M.) = 3(y~,-1) (2.039 


and the density function becomes 


w= pee 7) 

It 1s now possible to aoiee for the particle velocity, $6, 
by substituting equation (2.12) into the one dimensional con- 
tinuity equation and solving for the constant with $=$,, at 
—F=1 (the shock front location). The continuity equation is 


given by 


CAD iG Ry 
a eee BM SM 


(o-€) 
where £8 = R.R./R2 is the shock decay coefficient. With the 
above substitution and integration, the particle velocity 


profile is found to be 


@?= 9,6 


14 





By substituting the density and particle velocity pro- 


files into the momentum equation 


_£) 2% 1 9f ~ _gy 3d 
(OE) oe + BO + Foe = —BMea 


and integrating, the pressure profile is found to be 


y > a 
ff, + qa tb) 1) Ugh 8 (OMe at) 


where the constant of integration was solved with f=f, at 
g=l. : 
The shock. decay coefficient, 8, can be found using the 
assumption that shock radius is given by 
R,. = ent" (2.14) 
Experimental work done by NASA [Ref. 4] has found that 

N= 0.8 
for water. Substituting equation (2.14) into the expression 
for 8, the decay coefficient becomes 

Ge e-=0, 25 


It is now possible to solve the pressure distribution 


function given the shock Mach number, M.- 


C. SHOCK MACH NUMBER 

By considering the initial energy being released at the 
center of symmetry in the fluid it is possible to equate the 
initial energy to the integral energy density behind the 
shock. The energy liberated initially, E,, is distributed 


| such that 


R 
Eg = f (e-eg+3) o4nr*dr (2.15) 


es) 








Using the assumptions of a strong hemispherical shock as well 
as the density profile and particle velocity profile as de- 
rived previously, the integral equation (2.15) in non-dimen- 
Sional form becomes 


i 
E. = 2p .C2R2M2 fo c 


O O'S S42 aut 11) O  FOTE TY ED “as {2 16) 


The values of Por and n can be found for a given fluid. 


oO 
By choosing a Mach number the dependent variables $9, and jy, 
are known. Equation (2.16) can then be integrated and 
2E,/Rs versus shock Mach number can be plotted. Given the 
energy released, Eor it is then possible to find the shock 
Mach number as a function of shock radius from the plot of 
2En/R- versus M.. 

There 1s however, an anomaly at a shock radius of zero. 
To avoid a finite amount of energy in zero volume or an in- 


finite energy density, reference 5 suggests the initial 


radius be given as 


E 
R, = (—— 1/3 (2 19) 
PCG 27 


Equation (2.17) was then used as the starting point for the 


solution of the theoretical pressures behind the shock. 
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III. EXPERIMENTAL PROCEDURE 


The basic components of the ballistic range used for this 
study are shown in figure 1. A down-range view of the bal- 
listic range in figure 2 depicts the chronograph screen ar- 
rangement and test tank. The two spark sources above the 
tank were used as light sources for the shadowgraphs. Two 
rifles were used, a .22 caliber and a Remington .222 caliber. 
Commercial Remington .22 caliber Long Rifle ammunition was 
used. The .222 ammunition was hand loaded for testing so 
that the velocities of the projectiles could be adjusted and 
would be known more precisely. Sierra bullets classified as 
Spitzers were used. Figure 3 displays the corresponding 
Shapes of the projectiles. The mass of the .222 caliber and 
-22 caliber projectiles are listed in figure 4 along with 
the parameters used to calculate the energy of each configu- 
ration. The area and drag coefficient parameters listed in 
figure 4 were experimental values used to calculate the ve- 
locity decay that was experienced during projectile passage 
through the test tank. 

AVTRON, No. A914T333, chronograph screens were used to 
provide start and stop pulses to Monsanto 101 B counters. 

As the projectile passed through the screen the pulse circuit 
was energized. The first screen started the first counter 
and the second screen stopped the first counter and simultan- 
eously initiated the second counter. The third screen 


stopped the second counter. The projectile impact velocity 


waS computed using these times and distance measurements. 
i 
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Figure 4. Projectile Parameters 
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The events within the fuel cell were timed with two addi- 
tional counters. A fourth chronograph screen provided a 
Signal pulse which was distributed to the adjustable time 
delay unit (ATDU) and the two remaining counters to melee te 
them. The main element of the adjustable time delay unit 
circuitry was a Signetic, N74123, retriggerable, monostable 
multivibrator as shown in figure 5. Within the ATDU were 
two duplicate delay circuits to provide a 100 volt output 
pulse to each spark source. The time delay for both circuits 
in the unit was adjustable; however, the duplicate ae ore 
delay was initiated by the first circuit's output. This re- 
sulted in the delay time for the second output pulse from 
the ATDU being the sum of the delay time of the first cir- 
cuit plus the delay of the second circuit. The delay time 
range for the first circuit was 500 to 1600 micro-seconds, 
while the second circuit provided an additional delay of 5 
to 15 micro-seconds. The 100 volt output pulses going to 
the spark sources were also used to stop the Provecunter st 
The first output pulse went to the spark source providing 
illumination for a shadowgraph in the horizontal plane, 
while the second pulse triggered the spark source for a 
shadowgraph in the vertical plane. Each counter therefore 
measured the time the shadowgraphs were taken. 

Two 16-inch collimating mirrors were used to pass light 
through the tank. The spark sources were located at the 
focal point of their respective mirrors. Figure 6 demon- 


Strates the arrangement of the tank with respect to the 
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Figure 6. Schematic of Shadowgraph Apparatus 
(Down-range View) 
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Spark sources, mirrors and photographic plates. The side 
and bottom walls of the tank were made of plexiglass through 
which very little light attenuation occurred. One side wall 
and the bottom wall had grid lines etched in the plexiglass 
for use in measuring the shock radius and cavity distance. 
High speed Polaroid film, Type 57, was used for the shadow- 
graphs. A ground, first surface mirror was used to deflect 
the light in the vertical plane 90 degrees. 

The tank or test cell had front wall inner dimensions of 
17.5 inches wide and 17.0 inches high with a depth of 17.0 
inches. Figure 7 contains a view of the tank. The main mem- 
ber of the front wall was one-inch thick aluminum and served 
aS a mounting plate for the specimens used to simulate fuel 
cell wall material. To accommodate the test material the 
mounting plate had a circular cut-out 15.0 inches in diameter. 
The specimen was mounted behind the mounting plate with 16 
bolts arranged in a circular pattern. A front wall of 7075- 
T6 Aluminum, .090 inches thick, with a prepunched entry hole 
is shown in figure 8. This wall was used in the experiments 
with the exception of the Few tests which used a solid front 
wall. The rear wall of the tank was identical to the front 
wall. 

To provide support for the front and rear walls aluminum 
angle strips were placed longitudinally at the four corners 
of the tank. The angle strips were bolted to the front and 


“Leamemounting plates. In addition to providing structural 
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Figure 8. Front Wall Test Specimen, 
Pre-punched Hole 
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rigidity, the angle strips retained the one-inch thick plexi- 
glass tank walls. The top of the tank was open. 

Room Temperature Vulcanizing (RTV) rubber sealant was 
used at all joints to make the tank water tight. At the two 
highest energy levels the sealant around the front wall had 
to be replaced several times. Apparently the shock pres- 
sures were sufficient to deteriorate the sealant bond to the 
plexiglass and aluminum. Tape was used to prevent water 
spillage through the prepunched hole and to avoid any signif- 
icant energy loss during penetration. 
| During initial testing the tank was set on the top of 
the stand. It was discovered that although the tank was re- 
strained by the stand in longitudinal and lateral motion, 
angular tank movement was created by the moment arm and the 
force of impact at the center of the tank. This invalidated = 
the shock radius data and the tank was subsequently bolted 
to the frame. It was also discovered that the tank and 
stand would move approximately eight inches during a test 
firing. To alleviate this problem the front two legs of the 
stand were bolted to the floor. Eventually these bolts 
worked loose and finally all four legs were bolted to the 
floor and 150 pounds of sandbags were placed on the stand. 

The experiment was run at three energy levels. One en- 
ergy level was tested before continuing with the next higher 
energy level. By varying the time delay with the ATDU, 
shadowgraphs of the shock radius versus time were taken. It 


was also possible to measure the projectile distance into 
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the tank versus time. Figure 9 shows a pair of typical 
shadowgraphs. The impact velocity of the projectile was 
computed from the times on the first two counters and the 
distance between chronograph screens. The shadowgraph times 
were measured by the third and fourth counters. 

Acquiring data for the .22 caliber was found to be very 
difficult because of the deviation in impact velocities from 
One test firing to another. A good estimation of velocity 
was required prior to a test to compute projectile time of 
flight from the fourth chronograph screen to the tank wall. 
That time plus the desired time for photographing the shock 
radius determined the delay time used to set the adjustable 
time delay circuit. Since the .222 caliber projectiles were 
hand loaded, the impact velocities were more predictable. 
However, in the .222 caliber tests the velocity had a suffi- 
client deviation to make it difficult to shadowgraph a shock 
radius smaller than one inch. 

The data were reduced by first calculating the projectile 
impact velocity. The time of flight of the projectile from 
the fourth chronograph screen to the tank was then calcu- 
lated. This time was subtracted from the elapsed time of the 
shadowgraphs as measured by the third and fourth counters. 
The result was the time corresponding to the shock position 
from initial impact at the £rene waPeemonOckhmeragdius, meas— 
ured on the shadowgraph, was then plotted at the correspond- 
ing time. Measurements of the shock radius were accurate to 


HOoOeetneness 
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Figure 9. Shadowgraphs After Projectile Impact Through 
a Pre-punched Front Wall 





Projectile distance into the tank versus time data were 
also accumulated. The shadowgraph at 24 micro-seconds in 
figure 9 shows that the projectile was not distinctly visible. 
Work done by McMillen [Ref. 6] implies that the projectile 
nose is at the tip of the light cusp seen at the leading edge 
of the cavity shadow. Light refraction caused by density 
gradients in the flow field forms this nearly triangular re- 


gion called the cusp. 
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IV. RESULTS AND CONCLUSIONS 


This research was conducted to werity the theoretical 
prediction of shock radius as a function of time. In all 
but four cases, the data collected utilized projectiles im- 
pacting the fluid through a pre-punched hole in the front 
wall of the tank. There were four data points taken where 
the projectile passed through a solid front wall before com- 
ing in contact with the fluid. This was done at the conclu- 
Sion of the experiment to see if there was a dapamture from 
the trends of the previous data and to obtain a qualitative 
estimate of the extent of shock phase damage. The solid 
front wall data were taken using the highest energy level 


tested (12323 in-lb). 


A. SHOCK RADIUS VERSUS TIME - PRE-PUNCHED FRONT WALL 
Figures 10, 11, and 12 show the experimental data and 
the theoretical prediction of shock radius versus time. 
Seldom were projectile velocities prior to impact constant 
during the testing. Changes in velocity have a direct in- 
fluence on the initial energy. The data were plotted and 
compared, however, assuming constant impact energies. It 
was therefore necessary to determine the shock radius sensi- 
tivity to changes in energy. A plot, using theoretical pre- 
dictions of the shock radius versus impact energy for a 
CONS rane Linc LI gure I3jmmindicates that there as an insig- 


nificant change in shock radius for typical experimental 
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variations of impact energy about the reference energy. The 
impact energy calculated from experimental data of the indi- 
vidual tests was within +200 in-lb of the comparison energy 
level. Corresponding changes in shock radius were of an 
order of magnitude smaller than the accuracy of the radius 
measurements. Therefore, no corrections were applied to the 
measured shock radii for experimental differences in impact 
energy about the comparison energy level. 

The data for an initial impact energy of 1600 in-lb had 
the best agreement with theory. Theory assumes that the 
total projectile energy is deposited at the point of impact. 
The shadowgraphs indicate, however, that none of the projec- 
tiles came to rest at the point of impact. At the higher 
energy levels (7493 in-lb and 12323 in-lb) the data fell be- 
low the theoretical results which implies generation of 
weaker shocks than those predicted by theory. This may be 
explained by considering the projectile as a piston in the 
early stages of shock formation. It would be expected that 
the rate of change of fluid volume displaced by the piston 
determines the momentum exchange between the projectile and 
fluid, hence the shock strength. The shape of the .22 cali- 
ber projectile is noticeably more blunt than the .222 cali- 
ber projectiles (see figure 2). This means that the .22 
caliber projectile would have a greater rate of change of 
fluid volume displaced in the formation stage of the shock 
phase. This would produce conditions which are more closely 


represented by the assumptions of the theory. It is expected 
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that blunt projectiles or warhead fragments produce shocks 
that agree closely with the theory. 

Eventually the shock becomes acoustic because of geomet- 
rical attenuation. The theoretical predictions for the time 
at which the shock is acoustic were 22 micro-seconds at 1.71 
inches for the energy level of 1600 in-1b, 36 micro-seconds 
at 2.81 inches for the energy level of 7493 in-1lb and 44 
micro-seconds at 5.40 inches for the energy level of 12323 
in-lb. Theory assumed no counterpressure (P,<<P). For this 
reason it 1S expected that the shock will actually become 
acoustic earlier than theory predicts. The experimental re- 
sults for this study exhibited too much scatter to accurately 
determine when the shock becomes acoustic. It is reasonable 
to assume the duration of the shock phase is accurate in 
order of magnitude. The cavity begins to oscillate at a time 
on the order of 1000 micro-seconds, which implies a distinct 
time separation between the shock and cavity phases. This 
research has shown that the shock phase and cavity phase may 


be studied independently. 


B. SHOCK RADIUS VERSUS TIME - SOLID FRONT WALL 

The four data points taken for the case of a projectile 
passing through a solid front wall and then into the fluid 
are shown in figure 12 as X symbols. Time limitations pre- 
vented the collection of more data. The limited results do 
not indicate a significant departure from the pre-punched 
data. Figures 14, 15 and 16 show typical damage to the .043 


incghethuck 6061—-T6 Aluminum front wall specimen. The*front 


Sry 








Figure 14. Solid Front Wall After Test (6061-T6 Aluminum, 
Wosminemeoick, E> 12323 in-1b) 





Figure 15. Typical Damage to Solid Front Wall at Point of 
Impact (Front view, 6061-T6 Aluminum, .043 
thick, E, = 12323 in-1b) 
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wall experienced plastic deformation, became concave isuiie 
outward direction and buckled. There was some additional 
deformation at a few of the restraining bolt locations be- 
lieved to have been caused by a slightly unsymmetrical load 
distribution. This was a result of projectile impact not 
being at the true center of the front wall. Projectile pene- 
tration caused no visible signs of cracking; although there 
was permanent deformation, the pressure load on the front 
wall was not sufficiently severe to produce cracking. 

Soper [Ref. 7] calculated the percent energy change ver- 
sus impact velocity for penetration of various thickness alu- 
Minum plates with no fluid foundation. For the specimen and 
‘impact velocity used in this study there was just under a 10% 
loss of energy predicted in penetrating the front wall. Using 
the theoretical prediction of shock radius versus impact en- 
ergy in figure 13, there would be an insignificant change in 
the shock radius for such a 102 change of impact energy. 

The shape of the projectiles after penetrating a solid 
front wall was significantly different than the projectiles 
which passed through a pre-punched hole. This was the re- 
sult of the copper jacket on the projectile being deformed 
by both penetration and the stagnation pressure loading on 


the projectile. 


C. PROJECTILE vt DOC tah Gey nee lth er LU eo 
The impetus for investigating the projectile velocity de- 
cay in the fluid is the requirement to know the time rate of 


change of energy exchange in order to study the cavity phase 
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[Ref. lJ]. It is difficult to predict velocity decay because 
of the many dependent variables of the problem. The analy- 
sis of velocity decay suggested by Lundstrom [Ref. 1] was 


utilized. The equation of motion used was 


av 1 2 
D= -m— = —0,V°C.A 4.1 
ates 2” Sg D in 


Assuming a constant deformed projectile area, constant drag 
coefficient and constant projectile mass it was possible to 
rearrange equation (4.1) and integrate, which yields 


(i G- / 2m) x 
ys ee? (4.2) 


a Ae 


Oo 


Further rearranging and integration from initial conditions 
of x=0 at t=0 yields the following expression for projectile 


position as a function of time 


=). 
x = = £n(BV,t+1) (4.3) 
where 
P¢CpA 
B= —-—— 
2m 


The density for the and ieetebed fluid, H,O, was assumed con- 
stant. In calculating the theoretical results of equation 
(4.3) the value used for V, was the projectile velocity im- 
mediately prior to impact with the fluid. It may be more 
Suitable to use the velocity of the projectile at the comple- 
tion of shock generation phase as the initial velocity in the 
etbuid, Vo° Figure 17 is provided to point: out a possible 


discontinuity of the projectile energy and velocity profile 
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during shock formation. The initial location and time of 
this event would still be adequately described by x=0 and 
t=0. A value for the change in projectile energy due to 
shock formation, AE .s or the change in projectile velocity 
due to shock formation, AV. would make it possible to solve 
iene Vo: If the experimental data had less scatter than in 
this study it would have been possible to determine the value 
of AE, required for the theoretical results to be coincident 
with such experimental results. It would then be possible 

to solve for V,. 

Typical projectile shapes after coming to rest in the 
fluid are shown in figure 18 for each energy level tested. 
Figure 4 shows the values for the average final deformed 
area, drag coefficient and mass used to calculate the theo- 
retical curves in figures 19, 20 and 21. Holm [{Ref. 8] has 
shown that the projectiles possess dynamic stability and do 
not tumble, which contributes to the deformation experienced. 
Stagnation pressures on the order of 10'psi were found to 
exist on the projectile surface. These pressures are large 
enough to cause plastic deformation of the lead .22 caliber 
projectiles and the copper jacketed, .222 caliber projectiles. 

Approximately fifty-one percent (51%) of the .22 caliber 
projectile impact velocity was lost within 1.5 inches of ae 
front wall. The .222 caliber projectiles with an initial 
impact energy of 7493 in-lb only lost 17% of their impact ve- 
locity at 1.5 inches. Projectiles with an energy 12323 in-lb 
lost 25% at the same penetration depth. Percentage differ- 


ences of the energy lost at the same distance from impact 


43 








Figure 18. 





Typical Deformation of Projectiles After 

Coming to Rest in the Fluid Filled Tank - 
Pre-punched Hole Results (L. to R.: .22 

caliber, Eo = 1600 in-lb; .222 caliber, 

EL = 7493 in-1b; .222 caliber, E, = 12323 
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may be partially attributed to the rate of deformation of 

the projectiles. It is reasonable that the lead .22 caliber 
projectile deforms more rapidly than the .222 caliber jac- 
keted projectile. It appears that the highest energy (12323 
in-lb) projectile has reached a limit at which the projectile 
jacket and lead core cannot sustain the enormous stagnation 
pressure. No longer was the deformation similar to the two 
lower energy projectiles (figure 18). The drag coefficient, 
deformed projectile area and percentage of projectile veloc- 
ity decay at a fixed distance for the two lower energy levels 
appears consistent. At the highest energy level, values for 
these variables do not appear consistent. When there is ex- 
tensive deformation, as in the case of the highest energy 
level projectile, the assumptions appear to be overly simpli- 
fied. 

The concept of an equivalent flat plate area, E=CA, may 
be more advantageous than considering the drag coefficient 
and deformed projectile area separately. Also included in 
figure 4 is the value for the equivalent flat plate area, f, 
which best matches theory to experiment. This parameter re- 
duces the number of dependent variables by one and describes 
the projectile drag characteristics with one number. 

In reference 7, minimum impact velocities to penetrate 
various thickness plates are given. The experimental veloc- 
ity decay data were extrapolated by theory to predict pro- 
jectile velocities at the rear wall of the tank. These were 


compared with the minimum impact velocity to penetrate the 
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rear wall. After traveling the depth of the tank, 17 inches, 
the .22 caliber projectile (E, = 1600 in-lb) was predicted 

er nave a velocity of 0.394 ft/sec. Experimentally, the .22 
caliber projectiles were found approximately six inches from 
the front wall where equation (4.2) predicts a velocity of 

71 ft/sec. The .222 caliber projectiles with an energy of 
7493 in-lb were predicted to have a velocity at the rear wall 
of 290 ft/sec, which was 32% of the minimum impact velocity 
required for penetration. The experiments with the .222 cal- 
iber projectiles at an energy of 7493 in-lb were the only 
test cases which made a noticeable dent but no penetration 

of the rear wall. The highest energy level projectiles were 
predicted to have a velocity of 105 ft/sec and there were no 
Signs of projectile impact with the rear wall. These results 
seem to indicate that the exponential velocity decay is rea- 


sonable where projectile deformation is not extensive. 


D. PRESSURE VERSUS RADIUS AND TIME 

Figures 22, 23, 24 and 25 show theoretical predictions 
of the shock phase pressure pulse and the pressure profiles 
behind the shock at various times. Figure 22 is provided to 
demonstrate the order of magnitude of the shock pressure from 
initial impact until the shock becomes acoustic for an energy 
level of 1600 in-lb. The remainder of the pressure versus 
radius curves were limited to a range of 5 micro-seconds to 
30 micro-seconds. This information was used to estimate 
pressure conditions which will be encountered when further 


experimental studies are undertaken to verify the theory. 
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The pressure profiles behind the shock were not plotted at 
negative pressures although negative pressures behind the 
shock front were predicted in the theoretical results. Since 
water is nearly incompressible, cavitation would occur as the 
pressures became negative. Cavitation has been observed on 
film records of similar hydraulic ram studies done at the 
Naval Weapons Test Center, China Lake and provided for the 
author by Lundstrom. 

Immersion of pressure transducers in the fluid is being 
considered as a means of meaSuring pressures. Figures 26, 

27 and 28 were plotted to estimate the magnitude of the pres- 
sures experienced at various distances from the front wall 
and the duration of the pressure pulse. The plots indicate 
the need for a transducer with a very fast response time. 

It is expected that actual pressures will be somewhat 
lower than predicted, with the possible exception of blunt 
nose projectile data. There are two reasons for predicting 
lower pressure profiles. As explained earlier, ‘the shock 
radius versus time data implied a weaker shock than the theo- 
retical predictions. This was the result of the assumption 
in the theory that total impact energy was deposited in the 
fluid during shock formation. A weaker shock would result 
in correspondingly weaker shock pressures throughout the 
shock phase. The second reason pertains to the pressure pro- 
file as the shock Mach number approaches one. At lower shock 
Mach numbers, the assumptions of undisturbed fluid pressure 


in front of the shock being much less than at the shock 
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(P(<<P) will become invalid. This assumption should cause 
the theoretical predictions to be higher than experimental 


results. 
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V. RECOMMENDATIONS 


The scatter of the experimental data measuring shock 
radius versus time could be reduced by placing the fourth 
chronograph screen closer to the tank. New counters with 
high frequency, temperature compensated oscillators to 
provide a better time base are also recommended. 

The ballistic range should be modified to provide the 
capability of firing higher energy projectiles. Shock 
radius versus time and projectile distance Beevete in 
the fluid versus time should be studied at these higher 
energy levels. 

Fluid pressure at and behind the shock should be com- 
pared with theory at various energy levels. 

This study used lead or copper jacketed, lead core 
projectiles and no tumbling was observed. Further stud- 
ies using steel jacketed projectiles should be completed 
to test for deviations from the data of this study. This 
information would provide the fuel cell designer with 
data typical of ballistic threats encountered in combat. 

Strain gauges should be placed on the front wall spec- 
imens to determine the stresses incurred during the shock 
phase. Strain data could then be compared with SATANS 
computer code. 

Higher energy projectiles which produce cracking or 


catastrophic failure of the front wall specimen should 


ae 





be tested. These data would provide the designer with 
the critical energy level at which shock phase damage 


may result in an aircraft loss. 
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